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ABSTRACT 

We present the first broadband A = 1mm spectrum toward the z = 2.56 Cloverleaf Quasar, obtained 
with Z-Spec, a 1-mm grating spectrograph on the 10.4-meter Caltech Submillimeter Observatory. The 
190-305 GHz observation band corresponds to rest-frame 272 to 444 fim, and we measure the dust 
continuum as well as all four transitions of carbon monoxide (CO) lying in this range. The power-law 
dust emission, F v = 14mJy (i//240GHz) ' is consistent with the published continuum measurements. 
The CO J = 6^5, J = 8^7, and J = 9^8 measurements are the first, and now provide the highest- J 
CO information in this source. Our measured CO intensities are very close to the previously-published 
interferometric measurements of J = 7— >6, and we use all available transitions and our 13 CO upper 
limits to constrain the physical conditions in the Cloverleaf molecular gas disk. We hnd a large mass 
(2-50xl0 9 M Q ) of highly-excited gas with thermal pressure nT > 10 6 Kcm -3 . The ratio of the total 
CO cooling to the far-lR dust emission exceeds that in the local dusty galaxies, and we investigate the 
potential heating sources for this bulk of warm molecular gas. We conclude that both UV photons and 
X-rays likely contribute, and discuss implications for a top-heavy stellar initial mass function arising 
in the X-ray- irradiated starburst. Finally we present tentative identifications of other species in the 
spectrum, including a possible detection of the H2O 2o,2— ! Ti,i transition at A ros t = 303 /im. 
Subject headings: ISM: clouds, galaxies: ISM, stars: luminosity function, mass function, instrumenta- 
tion: spectrographs, techniques: spectroscopic 



1. INTRODUCTION 



Since the first measurements of CO in 

the p owerful IRAS gala xy FSC10214 at 

z=2.3 (jBrown fc Vanden Boutl j 19921) nearly two 
decades ago, the field of molecular line measurements in 
ear ly-universe galaxies has g r own steadily (see review 
bv ISolomon fc Vanden Boutl (|2005l )). Several tens of 
sources are now known, drawn from a wide variety of 
parent populations including IRAS galaxies, optically- 
selected quasars, submillim eter galaxies , and now 
Spitzer-selected galaxies (e.g. iFraver et al.l ([2008)). As 
with the submillimeter and millimeter-wave continuum 
measurements, observations of high-z molecular lines 
benefit from a negative K-correction: the spectral lines 
generally carry more power at higher frequency. In 
particular for carbon monoxide, the run of the CO line 
luminosity with J typically increases up to the mid-J 
transitions (J~5-10) for actively starforming galaxies. 
This CO rotational spectrum when measured in energy 
units thus peaks in the 200-500 /im regime, making 
millimeter-wave spectroscopy a powerful tool for study- 
ing molecular gas in the early universe. Measurements 
of the CO spectrum across its peak constrains the 
temperature, density and total mass of molecular gas. 
While observations of mid-J CO transitions in the local 
universe are hampered by poor transmission in narrow 
atmospheric windows, a millimeter- wave systems can 
have approximately uniform sensitivity across the peak 
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of the CO SED for z >1. 

All high-z millimeter-wave observations to date have 
employed heterodyne receivers, originally with single- 
dish telescopes but soon after with large interferometer 
arrays. SIS receiver sensitivities, T rec ~ 30 K, are close 
to the photon background limit at mountaintop sites, and 
are an excellent match to the large-collecting area arrays. 
However, their instantaneous fractional bandwidths are 
at most ~ 10%, even with the new wideband backends 
(e.g. 8 GHz for ALMA). They thus require a separate 
tuning for each spectral line, and searching for unknown 
lines and/or unknown redshifts is time-consuming. A 
more comprehensive approach to finding redshifts and 
probing molecular gas contents of distant galaxies can 
be obtained with complete spectral coverage, to the ex- 
tent possible given the telluric windows. CO is but one 
of the important coolants, and it is generally accepted 
that other species, notably water, become important as 
the molecular gas becomes more excited. Z-Spec is the 
first spectrograph designed specifically for this type of 
measurement at wavelengths longward of the mid-IR; a 
first-order grating with simultaneous coverage over the 
entire 1 mm atmospheric window. 

After a brief introduction to the Cloverleaf quasar, Sec- 
tion [2] provides a description of Z-Spec and our observa- 
tions at the CSO. The line and continuum parameters 
extracted from our spectrum are presented in Section [3j 
and the findings are discussed in Section |4] 

1.1. The Cloverleaf 

We have observed the lensed z = 2.56 broad- 
absorption-linc quasar system H1413+1143, also known 
as the Cloverleaf. The source was originally detect ed in 
an optical spectroscopic survey (Ha zard et al.lll984l ). but 
a warm dust component also emits some ~ 7 x 10 13 L© 
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Fig. 1. — Z-Spec Spectrum toward the z=2.56 Cloverleaf Quasar. Only one channel of the total 160 has been removed due to high 
noise, at 209.0 GHz — a value near the baseline is plotted in its place, noted with a diamond. A fit to the continuum and six spectral lines 
is overlaid with a light red curve. Red squares with error bars and spectral widths denote continuum measurements from MAMBO and 
Plateau de Bure (Section 14.11 1. Thin dot-dash lines mark transitions frequencies for which upper limits are extracted . The J = 7 — > 6, 



J = 8 -» 7, and J = 



CO frequencies are marked in green with 13,7, 13,8, and 13,9. Frequencies of the two water transitions with 



Eupper < 140 K in the band are marked in cyan. The LiH and CH + frequencies and a tentative absorption feature are are marked with 
magenta. Upper limits are provided for all transitions except 1 3CO J = 9^8 in Table [T] See also section [1.5.21 

Subsequent observations of the J = 3 — > 2 and J = 7 — > 6 
transitions with the BIMA, OVRO and IR AM inter- 
ferometers have steadily improved t W ilner et al.l 119951 : 
lYun et al.lll997D . culminating in an IRAM J = 7 — > 6 
map which fully r e solves the four Clover leaf compo- 
nents |A1 oin et al] (Il997f) [her e after A97], iKneib et all 
(|1998al) ). IVenturini & Solomon! (|2003l ) used the J = 7^6 
image to model the CO J = 7 — > 6 source and lens and 
find (after correcting to the modern cosmology) that the 
intrinsic source is a disk with radius of 650 pc, inclined by 
30°. The lens is formed with two galaxies of comparable 
mass at 0.25" and 0.71" from the line of sight. Our Z- 
Spec measurement complements this spatial information 
by providing a survey of several CO transitions spanning 
the peak of the spectrum. The uniform calibration allows 
us to anchor the total molecular gas conditions and ener- 
getics. The result provides constraints on the gas heating 
sources: stars or the active nucleus, and on the impact 
of the active nucleus on the putative host starburst. 

2. Z-SPEC INSTRUMENT AND OBSERVATIONS 



(intrinsic) in the rest-frame mid-IR ijLutz et al] [2007). 
ascribed to reradiated energy from the accretion zone. 
The SED also suggests a second distinct dust emission 
component dominating the rest-frame far-IR and sub- 
millimeter, but with a total bolometric luminosity ~5- 
10% that of the AGN dust component (jBarvainis et al.l 
H99l[l995T : IWei6 et al.ll2005h . The measurement of PAH 
features with intensities consistent with a starburst in 
this far-IR / submm component leads to the picture that 
the Cloverleaf is a composite object: an usually powerful 
partially-embedded AGN accompanied by a weaker but 
still tremendous starbust which alone wou ld be similar 
to a submillimcter galaxy (jLutz et al.ll2007l ) . These large 
luminosities, combined with a substantial lensing magni- 
fication factor (~ 11) makes the Cloverleaf an excellent 
laboratory for studying the AGN / host interaction in 
what is likely the era of peak activity in galaxies. In par- 
ticular, the molecular gas reservoir is the raw material for 
both the star formation and ultimately the nuclear ac- 
cretion, and the impact of UV and X-ray photons from 
the stars or AGN on the molecular gas is a key aspect of 
the starburst / AGN interaction. 

Some of the earliest high-z molecular gas measure- 
ments were detections of th e CO J = 3— >2 transition red - 
shifted to the 3 mm band (|Barvainis et al.lll994 I1997D . 
These measurements confirmed the presence of a large 
gas reservoir (~ 10 10 M Q ) inferred from the dust SED. 



Z-Spec is a single-beam grating spectrometer which 
disperses the 190-308 GHz band across an linear array of 
160 bolometers. The grating approach is novel: a curved 
grating operates in a parallel plate waveguide which is fed 
by a single-mode corrugated feedhorn. The instrument 
coupling to the CSO telescope is thus well-approximated 
with a Gaussian-beam approach, and the measured ef- 
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ficiencies and beam sizes are consistent with this. De- 
tails of the grating design and testing can be found in 
iNavlor et alT(l2003T); iBradford etafl ((200l ; lEarle et all 
(|2006fk llnami et all (|2008|) Thedetector spacing in- 
creases from ~ 1 part in 400 at the low-frequency end 
(Av = 500 MHz), to ~ 1 part in 250 at the high- 
frequency end of the band (Av = 1200 MHz), while the 
spectrometer resolving power runs from ~ 1 part in 300 
at low frequencies to ~ 1 part in 250 at the high fre- 
quency end of the band. Thus the system is marginally 
under-sampled, especially at the high-frequencies. 

Spectral profiles for all 160 channels have been mea- 
sured with a long-path Fourier-transform spectrometer, 
with channel center frequencies adjusted slightly per ob- 
servations of multiple transitions in the spectral standard 
IRC+10216. Refinement of the Z-Spec frequency scale is 
ongoing as we incorporate line measurements from an 
increasing library of astronomical sources; for the data 
presented here, we are confident that the channel fre- 
quencies are known to better than 200 MHz across the 
band, or approximately one fifth of a channel width. We 
expect to improve this further in the near future. 

The entire structure is cooled with an adiabatic de- 
magnetization refrigerator and operates at temperatures 
between 60 and 85 mK to facilitate photon-background- 
limited detection. The bolometers developed at the 
JPL Microdevices Laboratory are individually-mounted 
silicon-nitride micro-mesh absorbers with quarter-wave 
backshorts, read-out with neutron-transmutation-doped 
Germanium thermistors. With an operational opti- 
cal loading is ~ 1 — 3 x 10~ 13 W, and phonon NEPs 
of 4 x 10~ 18 WHz -1 / 2 , these detectors are the most- 
sensitive, lowest-background bolometers fielded to date 
for astrophysics. 

Z-Spec observes in a traditional chop-and-nod mode, 
with the secondary chopping between 1 and 2 Hz, and 
a nod period of 20 seconds. Because it is not possible 
to modulate the spectral response of the instrument rel- 
ative to the bolometer array and the spectral resolution 
elements are not oversampled, it is critical to both in- 
sure excellent array yield, and carefully calibrate each 
detector's response. The situation is complicated by the 
fact that the both the bolometer loading and the bath 
temperature are changing throughout a typical observ- 
ing night, while primary calibration sources are relatively 
scarce. To address this, we have built a library of plan- 
etary observations in varying conditions, and we fit the 
dependence of each bolometer's response on its operat- 
ing voltage, a proxy for the combination of bath tem- 
perature and optical loading. This relationship provides 
a calibration correction which is used to bridge intervals 
between astronomical calibration observations. Based on 
the self-consistency obtained with this scheme on planets 
and quasar calibrators, we estimate that the channel-to- 
channcl calibration uncertainties are less than 10%, ex- 
cept at the lowest frequencies which are degraded by the 
wing of the 186 GHz atmospheric water line. 

The Cloverleaf was observed on 2008 April 8 and April 
15, with a total observing time (including chopping) of 
7.9 hours, equally split between the two nights. On April 
8, the weather was excellent, with T225GHZ between 0.04 
and 0.06, while April 16 had higher opacity: T225GHz~ 
0.15-0.2. Each channel's signal is summed individually 



on a nod-by-nod basis, with each nod individually cali- 
brated and weighted by the inverse square of the detector 
noise (in Janskys) measured around the chop frequency. 
No modifications are made to the overall spectral shape 
at any point. Z-Spec's sensitivity is weather dependent: 
on April 8, the sensitivity per spectral channel was 300- 
500 mJy v / sec, but degraded to 500-700 mJy -/sec on the 
15th, with greater degradation at the band edges, where 
the opacity increase is larger. The measured sensitivi- 
ties are consistent with a simple model in which photon 
noise from the sky and telescope is the dominant contri- 
bution to the system noise. The final RMS uncertainty 
in each channel is shown with the error bars in Figure 1. 
It ranges from 2.5 to 4 mJy for most channels, with a 
median of 3.3 mJy. 

3. RESULTS, CONTINUUM AND LINE FLUX 
EXTRACTIONS 

The Z-Spec spectrum toward the Cloverleaf system is 
shown in Figure [TJ To extract line fluxes, we perform 
a simultaneous fit to a single component power-law con- 
tinuum and the multiple lines. Each channel's measured 
spectral response profile is used in the fitting since the 
spectrometer is not critically sampled, and the grating 
has spectral sidelobes at the ~1% level. All line widths 
are fixed, since Z-Spec is not sensitive to widths be- 
low ~1000 kms' 1 . The PdB measurements of the CO 
J = 3^2 sh o ws a w idth of 420 kms _1 (W03), while the 
iKneib et alj (|1998aD spectrum of J = 7 — > 6 has a Gaus- 
sian fit to some of the profile giving FWHM^450 km s~ 1 , 
though the profile is non-Gaussian and the true RMS 
width of the profile is larger than this. It is conceivable 
that the line width increases with J — this would be ex- 
pected if warmer material lies preferentially closer to the 
point-like nucleus. We adopt 500 kms" 1 for the Z-Spec 
fits. This choice is not crucial; we have found that the fits 
for the integrated line fluxes are not strongly sensitive to 
the adopted line width for values below ~800 kms -1 . A 
total of six spectral lines are fit, the results are presented 
in Table Q] 

The [CI] J = 2 -> 1 (y = 227.5 GHz) and the 
CO J = 7 — > 6 (v = 226.7 GHz) transitions are sepa- 
rated by only 1000 kms" 1 , or about one Z-Spec chan- 
nel, so blending is a problem. To address this, we adopt 
an iterative approach to make appropriate use of prior 
information: we first fit the continuum plus the three 
unblended CO lines (J = 6^5, J = 8^ 7, and J = 9 -^8) 
to determine their fluxes and a CO redshift. We find 
z = 2.5585 ±0.0015, a value in good agreement with pre- 
vious measurements (our error derives from the accuracy 
of our frequency calibration, and can be improved in the 
future). With the redshift determined, we then fit fluxes 
to the CO J = 7 — > 6 and [CI] J = 2 — > 1 lines, find- 
ing 45.3 ± 6.4 and 8.5 ± 7.4 Jy kms -1 , respectively (see 
Table [I]). Our current frequency scale uncertainty corre- 
sponds to interchanging ~4 Jy km s" 1 of flux between the 
two transitions about these fit values, with the sum of the 
two conserved, a value which is less than but comparable 
to the statistical errors. Our results are consistent with 
the various J = 7 — * 6 integrated intensity measurements 
from the Plateau de Bure. B97 and A97 quote 43.7 ± 2.2 
and 43.3 ± 2.4 Jykms -1 , respectively, using a fit to a 
Gaussian profile with AVfwhm = 376 kms -1 (adopted 
based on the lower- J measurements). However, A97 also 
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report 50.1±2.8 Jykms" 1 when the linewidth is fit as a 
free parameter (finding AVfwhm = 480 kms" 1 ). Given 
the linewidth discussion above, we adopt the larger of 
these as our input to the analysis, but adopt a 15% sys- 
tematic uncertainty to accommodate the range of mea- 
surements. As the CO temperature spectrum plotted 
in Figure [2] shows, this is in good agreement with the 
J = 6— >5 and J = 8—>7 measurements, giving confidence 
to the overall calibration and noise estimates. Our [CI] 
flux is also formally consistent with the W05 measure- 
ment (5.2 ±0.3 Jykms- 1 ). 

An emission feature at 277.6 GHz is also fit, tenta- 
tively identified as water and discussed in Section f4. 5. II 
Finally, we have obtained upper limits to six other tran- 
sitions, discussed briefly in Section [4.5.21 The potential 
emission feature at 208 GHz is unidentified at present. 

4. DISCUSSION 

4.1. Dust Continuum 

The best fit continuum across the Z-Spec band is 
F v = 14.1 ± 0.4 mJy(i7240GHz) 3 - 91±0 - 17 . This is con- 
sistent with previous measurements, for which the best 
compilation is in W03 (their Figure 3); we show the lat- 
est measurements in our band in Figure [T] The MAMBO 
bolometer measurement lies very close to the Z-Spec 
spectrum; the PdB continuum measurement at 224 GHz 
is slightly lower, but we note that this data point lies be- 
low the best fit of W03 including all the higher-frequency 
data. Indeed, the various measurements and limits from 
97 GHz to 400 GHz do not all lie on a single greybody 
curve, but our best-fit exponent of 3.91 from 190 to 300 
GHz is consistent with the W03 best fit: dust at 50 K 
with (3 = 2 and with unit emissivity between 50-300 /im. 

4.2. CO Excitation and Radiative Transfer Modeling 

The run of CO line luminosity with J for all pub- 
lished measurements is shown in Figure [2] The lumi- 
nosity scale shows the quoted statistical uncertainties for 
the various measurements. Systematic errors are clearly 
present judging from the published J = 3 — > 2 through 
J = 5 — > 4 measurements, and we have assumed that 
systematic effects produce a normally-distributed uncer- 
tainty of cr = 20% of the measured flux for the non-Z- 
Spec measurements. This statistical uncertainty is shown 
in the temperature scale. 

Our approach is to consider the properties of the gas 
which is confined to the physical size inferred by VS03, a 
650 pc-radius disk inclined by i = 30° to present a 559 pc 
semi-minor axis. VS03 resolve 78% of the total J = 7^6 
flux in this source, and we adopt this fraction for other 
mid- J transitions, in particular for J= 8—^ 7, our adopted 
reference transitions for bolometric calculations given its 
clean measurement with Z-Spec. The large observed line 
luminosities arising from the VS03 disk immediately im- 
ply that the molecular gas in the disk is very warm. If 
the disk radiates isotropically with a velocity FWHM of 
500 kms _1 ,the brightness temperatures must be Trj = 
37, 40, 31, 20 K, for the J = 6^5, J = 7^6, J = 8->7, 
and J = 9 — > 8 lines, respectively. The corresponding 
physical temperatures assuming optically-thick, thermal- 
ized emission are even higher: 50, 56, 49, 39 K. An area 
filling factor or optical depth less than unity would cor- 
respond to even higher physical temperatures. Adopting 
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Fig. 2. — Cloverleaf CO spectral line energy distribution (SLED) 
using all available transitions. The top panel is in bolometric en- 
ergy units, the bottom in brightness temperature units, referred 
to the VS03 source size assuming m = 11. We have adoped a 
20% systematic uncertainty for the measurements other than those 
from Z-Spec, shown in the temperature plot. A model correspond- 
ing to conditions near the peak in the pressure likelihood (T = 



56K,n H2 = 2.8 X 10 4 ,cm- 3 , N co /dv = 1. 



10 6 ( 



; km 



normalized to match the observed J = 8 — > 7 flux) is overplotted as 
a solid line in both scales. Dashed lines in the temperature scale 
show themalized blackbody emission for temperatures of 25, 50 
and 100 K (higher T is shallower slope), arbitrarily normalized to 
the model value of 39 K at J = 6— > 5. To produce 78% of the total 
observed intensities, the area filling factors relative to the VS03 
disk would have to be 3.3, 1.1, and 0.46, respectively. 

a velocity width smaller than the 500 kms -1 that we 
assume would also correspond to even higher physical 
temperatures. 

For a more quantitative a nalysis, we begin w i th a y ari- 
ant of the RADEX code (|van der Tak et all [2007h to 
model the CO excitation and radiative transfer. RADEX 
is a non-LTE code that uses an iterative, escape prob- 
ability formalism to treat the line radiative transfer. 
Three independent physical variables provide the in- 
put into RADEX: the gas density (nnj, the kinetic 
temperature (T), and the CO column density per unit 
linewidth Nqo/Av, which sets the optical depth scale. 
RADEX calculates the excitation temperatures, line op- 
tical depths, and line surface brightnesses in the CO lines. 
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TABLE 1 

Far-IR and Submillimeter Transitions Observed in the Cloverleaf 



transition 


Vrest [GHz] 


Euppcr [K] 


!/obs [GHz] 


flux [Jy km/s] 


err [Jy km/s] a 


L [10 s L ] b 


reP 


CO J = 3->2 


345.79 


33.2 


97.2 


13.2 


1.7 


0.55 


W03 


CO J = 3^2 


345.79 


33.2 


97.2 


9.9 


0.6 


0.41 


B97 


[CI] J = 1^0 


492.16 


23.6 


138.3 


3.6 


0.4 


0.21 


W05 


CO J = 4^3 


461.04 


55.3 


129.58 


21.1 


0.8 


1.17 


B97 


CO J = 5 — +4 


576.27 


83.0 


161.96 


24.0 


1.7 


1.66 


B97 


CO J = 6^5 


691.47 


116 


194.3 


37.0 


8.1 


3.1 


this work 


CO J = 7^6 


806.65 


155 


226.7 


50.1 


2.2 


1.9 


B97, A97 


CO J = 7^6 


806.65 


155 


226.7 


45.3 


6.3 


1.4 


this work 


rOTl 7 O ^ 1 

[U1J J — z — > 1 


ouy.o4 


37 


226.6 


5.2 


0.3 


U.oU 


W Uo 


[v^rj J — Z — r 1 




37 


226.6 


8.5 


7.4 


u.oo 


this work 


CO J = 8^7 


921.80 


199 


259.0 


51.4 


1.7 


5.7 


this work 


CO J = 9^8 


1036.91 


249 


291.4 


41.8 


5.8 


5.2 


this work 


far-IR dust 




115 K, peaking at 100 fim rest frame 




5.4xl0 4 


W03, L07 






Tentative Identifications and 2.5 <r Upper Limits 






H 2 2 ,2-li,i 


987.93 


101 


277.6 


20.3 


6.1 


2.4 




H 2 2i jl ->2 ,2 


752.03 


137 


211.4 


<14.2 


5.7 


<1.3 




Absorption feature 






288.3 


-17.4 


5.8 






LiH J = 2^1 


887.29 


21.3 


249.3 


<12.3 


1.9 


<1.9 




CH+ J = 1^0 


835.07 


43.1 


234.7 


<14.0 


5.6 


<1.6 




13 CO J = 7^6 


771.19 


148 


216.81 


<11.7 


1.7 


<1.1 




"CO J = 8-+7 


881.27 


190 


247.76 


<14.8 


5.9 


<1.6 





a Statistical uncertainties only, additional systematic uncertainties are included in the analysis. Luminosities are corrected for 
lensing, usin g a magnification fac tor of 11. c References — W03: IWeifi et al.l i2003|) , W05: IWeifi et aT]||2005H , B97: IBarvainis et all 
1119971) . A97: lAlloin et ail U997T) . 



Because of the uncertainties in the magnification and the 
beam filling factor (assumed to be the same for all tran- 
sitions), we compare the line ratios (with respect to the 
J = 8^7 line), rather than the absolute line fluxes. We 
have used RADEX in this way to generate the CO line in- 
tensities on a large grid in density (10 2 -10 8 cm -3 ), tem- 
perature (0-300 K), and CO column per unit linewidth 
(lO^-lO^cm^km-is). 

The CO column density per unit linewidth is not 
strongly constrained by the 12 CO line ratios alone, 
though very small optical depths are not allowed because 
some degree of radiative trapping is required to populate 
the levels up to J = 9 (T= 249 K). To address this, we 
also compare the upper limits for the 13 CO J = 7 — > 6 
and J = 8—>-7 transitions in our band with the output of 
a RADEX model for this species, assuming that its frac- 
tional abundance relative to 12 CO is 1/40 , a value de- 
rived in a multi- level study of the NGC 253 lHenkel et all 
(jl993h . 

We then follow the methodology of I Ward et "all (|2003l ) 
to generate likelihood distributions for nn 2 , T, and 
Nco/Av by comparing the RADEX results with the ob- 
served line ratios. We assume uniform priors in the log- 
arithm of all physical parameters, except for two prior 
constraints. First, the total mass of gas producing the 
CO lines in the disk cannot exceed the dynamical mass 
implied by the observed velocity spread and the size of 
the disk. This eliminates very large CO optical depths 
according to: 



A^CO MdynAcO 

dv 



1 



1 



1.4mr 



7ri?j cos i Av 



1.3xl0 18 cm -2 km -1 s. 

(!) 

Here the dynamical mass in the disk is based on the 
650 pc modeled size, the measured aspect ratio implying 
i = 30°, and a circular velocity of 375 km s -1 /sin 30°, per 
the B97 HWZI profile width (and confirmed by W03): 



M dyn = 8.5 x 1O 1O M . X co is the CO abundance rela- 
tive to H2 (taken at 2 x 10~ 4 ) and /i = 1.4 is the mean 
molecular weight in units of Mh 2 ■ 

The second constraint requires a minimum tempera- 
ture sufficient to produce the observed luminosity in the 
finite-sized disk, described above. Figure [3] shows the 
effect on the likelihood distributions of adding this con- 
straint: with it, the physical temperature is required to 
be above ~50 K, without it, the temperature likelihood 
suggests T~30-60 K. The constrained-temperature like- 
lihood suggests somewhat lower densities, and provides 
a better-defined total thermal pressure than the uncon- 
strained likelihood, peaking between 0.8-3x 10 6 K cm -3 . 
We note that the derived likelihoods depend on the as- 
sumed linewidth only through this second constraint. 
In its absence, the likelihoods depend only on the in- 
tegrated intensities. As described above, the adoption 
of 500 kms -1 linewidth value is actually a conservative 
value in this context: smaller linewidths corresponds to 
even higher temperatures, which would provide a more 
stringent constraint on the final likelihoods. 

While the disk dominates the mid-J line emission, it 
may not dominate the integrated large-beam low-J (e.g. 
J=l— >0 through J = 3— >2) fluxes. We explore the de- 
gree to which a single component model can account for 
all the transitions by both including and excluding the 
J = 3 — > 2 measurement (the lowest-J transition avail- 
able). The results are shown in Figure [4] Not surpris- 
ingly, higher excitation conditions are favored when the 
J = 3 — > 2 is excluded, but the effect is small, and in- 
dicates that a single component is a suitable model for 
all the observed transitions. Our adopted likelihoods for 
the physical conditions are thus shown with the black 
curves in both Figures [3] and [H yielding a temperature 
of 50-100 K and a density n H2 > 3 x 10 3 cm~ 3 . As an 
example, the spectrum in Figure [2] is overlaid with the 
fluxes predicted by a model near the peak of the pres- 



G 



Bradford et al. 





5 6 7 8 9 10 16 17 18 5 6 7 8 9 10 16 17 18 

Log P/k (cm- 3 K) Log N(C0)/AV (cm 2 km 1 s) Log P/k (cm- 3 K) Log N(C0)/AV (cm 3 km 1 s) 



Fig. 3. — Likelihood distributions for the physical conditions 
in the Cloverleaf host, with (red) and without (black) a prior ap- 
plied to insure that the temperature is sufficient to produce the 
observed luminosity in region with size given by the modeled disk. 
For both, the constraint in column density per linewidth is applied, 
eliminating Nco/At> values much above f0 18 cm~ 2 km _1 s. Oth- 
erwise, uniform priors are adopted in the logarithm of n, T and 
Nco / A«. All transitions are included. 

sure likelihood with T = 56 K, nn 2 = 2.8 x 10 4 cm -3 , and 
./Vco/Au = 1.6 x 10 16 cm~ 2 km -1 s. We note that for our 
peak likelihood conditions, the gas density and the opti- 
cal depth parameter are consistent with the gas having 
at least enough velocity dispersion to correspond to viri- 
alized motion under its own self-gravity . Rearranging 
Equation 2 in iPapadopoulos et alT (|2007l ) for our units 
gives: 
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where Nco / dv is in units of cm km s, Xco is the CO 
abundance relative to H 2 , a is a parameter ranging from 
1-2.5 per the cloud density profile, and the inequality ac- 
counts for the fact that the gas can of course be subject 
to more than its own gravity. For a density of 10 4 cm -3 , 
the inequality is satisfied as long as Nco /dv is less than 
3.8 x 10 18 cm -2 km -1 s. Our derived likelihoods thus im- 
ply velocity dispersions which exceed the virial require- 
ment by about an order of magnitude, perhaps due to 
large-scale turbulence and/or the influence of additional 
mass (of at most ~3 times that in the gas itself). 

4.3. Mass of Molecular Gas 

Our results are broadly consistent with the analysis 
of B97; cf. their Figure 2. As in that paper, we can 
directly compute (a likelihood for) the mass by scaling 
from the observed line luminosities and using the implicit 
information on the line optical depth from RADEX, in- 



FlG. 4. — Likelihood distributions for the physical conditions 
in the Cloverleaf host, both with (black) and without (red) the 
J = 3 — » 2 measurements included. The temperature and optical 
depth constraints are applied, but otherwise uniform priors are 
assumed in the logarithms of n, T and Nco/Aw. 

stead of relying on an uncertain CO X-factor The cal- 
culation depends somewhat on the cloud geometry, as it 
determines the ratio of mass to emitting area; we have 
assumed spherical clouds, for which the working relation- 
ship is 
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where N 17 = [N/ Av]/[10 17 cm~ 2 / (kms^ 1 )}, the inten- 
sity I is in cgs integrated intensity units (erg s _1 cm -2 
sr _1 ). This conversion results in a similar mass when 
applied to the various mid- J transitions, and Figure [5] 
shows the likelihood distribution for the total gas mass, 
again with and without J = 3 — >2. Our likelihood corre- 
sponds to a true molecular gas mass in the source of M ~ 
0.2 — 30 x 1O 1O M , consistent with previous gas mass esti- 
mates derived from CO, dust, and neutral carbon which 
range from 1.3-1. 6xl0 10 M Q (W05). We note for com- 
pleteness the consistency around the peak values in the 
mass and optical depth likelihoods. The gas column den- 
sity is the estimated M = 6 x 10 9 M Q distributed over the 
(projected) VS03 disk: N H ~ 4.6 x 10 23 c m - 2 . The CO 
optical depth parameter should be the corresponding CO 
column (N co = 2 x 10" 4 N H ~ 4.6 x 10 19 cm" 2 ) divided 
by the linewidth of the source (400-500 kms -1 ), giving 
N C o/dv- 0.9- 1.2 x 10 17 cm" 2 km -1 s. 

4.4. Origin of the Warm Molecular Gas 

The CO emission observed in the Cloverleaf is unusu- 
ally intense, even if the entire far-IR luminosity is pro- 
duced by a starburst with Lsb = £40- 1207ml ~ 5.5 x 
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Fig. 5. — LEFT: Likelihood for the total molecular gas mass 
traced with CO in the Cloverleaf host, assuming a magnification 
factor m of 11. RIGHT: Likelihood distribution for the total lumi- 
nosity in the all the CO transitions. Color coding is as in Figure [4] 



1O 12 L as suggested bv iLutz et all (|2007l ). W03. The 
most luminous CO lines, the J — 8^7, J = 9 — ^8 each 
emit a fraction ~ 10~ 4 of this total starburst luminosity 
(see Table [1]), and we measure the total CO luminos- 
ity through the likelihood formalism, finding that it is 
well-constrained at 4.6-6.8 x the J = 8 — > 7 luminosity 
(Figure [5]). We therefore have 

L co « 3.3 x 10 9 L a 6.1 x 1CT 4 L far _ m , (4) 

assuming that the same value of m applies to both the gas 
and dust. This is more extreme than for the nuclei of the 
nearby starburst galaxies. In NGC 253, for example, the 
CO SED (in energy units) in the central 180 pc appears 
to peak at J = 6— >5 or J = 7^ 6; and each of these lines 
carry 2-3 x 10 -5 times the far-IR luminos ity, so the total 
CO luminosity fraction is ~ 1-2x1 0~ 4 ([Bradford et al.l 
120031 iHailev-Dunsheath et alJl2008t iGiisten et alJl2006[ l . 
Moreover, rec ent mid-J CO obser vations with the ZEUS 
spectrometer ( Sta cev et al.l [2007T) in a sample of local- 
universe LIRGS and ULIRGS show mid-J (J = 6 — > 5 
and/or J = 7^ 6) CO fractional luminosities comparable 
to or less than those in NGC 253 (T. Nikola, personal 
communication). Thus molecular gas cooling relative to 
the far-IR dust emission in the Cloverleaf host demon- 
strably exceeds that in even the extreme local-universe 
systems. We now consider potential heating sources for 
this bulk of warm molecular material. 

4.4.1. Stellar Ultraviolet (PDRs) 

Rec ent PDR model calculations (Kaufma n et al.l 
(2006), hereafter K06) compute CO line intensities which 
can be compared directly with the observations. These 
models parametrize the space of physical conditions with 
two quantities: the incident far-UV flux (Go), and the 
molecular hydrogen density of the nascent molecular 
cloud (nn 2 ). They indicate that for relatively high den- 
sities (n > few x 10 4 ) and modest UV fields (1 < G?o < 
10 4 ), the CO J = 6 — > 5 and J = 7 — > 6 transitions can 
have intensities more than 10~ 4 times the total far-IR lu- 
minosity. We compare two measured ratios with model 
calculations available online as part of the K06 model: 
1) the measured intensity ratio of the CO J = 7 — > 6 to 
the total far-IR (= 8.8 x 10 -5 , assuming that they come 
from the same physical region), and 2) the CO J = 6^5 
to J = 2^1 (=23, assuming J = 2^1 is thcrmalized at 
the same temperature a sJ = 3^2). We find that these 



ratios can be consistent with the K06 PDR model with 
Go=l-5xl0 3 and uh 2 =l-4x 10 5 cm -3 . The density es- 
timate, when combined with our likelihood analysis of 
the thermal pressure, suggests gas temperatures in the 
range of 10-300 K, reasonable given the range of tem- 
peratures expected in the PDR. 

However, the neutral carbon intensities reported thus 
far are are not consistent with this picture. For the PDR 
parameters discussed above, the K06 model predicts that 
the [CI] J = 1 -> intensity relative to the CO J = 3 -> 2 
and J = 4— >3 intensities should be 0.2-0.4 and 0.1-0.3, 
respectively; ratios which are close to those observed. 
But the [CI] J = 2 ^ 1 to J = 1 ^ intensity ratio is pre- 
dicted to be ~6, while it is observed at only 2-2.5. The 
PDR model is predicting warmer neutral carbon than is 
observed, a fact which is generally difficult to reconcile 
with the conditions derived from the CO spectrum and 
which is discussed further in 14.4.21 below. 

For both total far-IR flux, and CO J = 7 -> 6, the 
observed intensity emerging from the disk is ~50 times 
greater than predicted by the single PDR model which 
fits the ratios. If UV photons are responsible for the heat- 
ing of the gas, then on average in the Cloverleaf disk, we 
are seeing ~50 PDR surfaces superposed along the line of 
sight. Correspondingly, the column density correspond- 
ing to a visual extinction (Ay) of 10 magnitudes is ~l/50 
of the total gas column estimated in Section 14. 2\ if it is 
averaged over the 650-pc disk. The derived far-UV field 
of Go ~l-5xl0 3 is modest, suggesting that the Clover- 
leaf PDR is not directly adjacent to an OB star as in an 
ionization-bounded HII region, but that the PDRs are 
typically 0.1 to 0.5 pc from the OB stars. This is quite 
reasonable, if the 5.5 x 10 12 L Q attributed to the star- 
burst arises from 100-10000 L© stars in a region of size 
given by the VS03 model, then the average interstellar 
separation is 0.2-1 pc. 

4.4.2. Hard X-Rays from the Active Nucleus 

The hard X-ray continuum of the active nucleus it- 
self could power the observed CO emission. Hard 
(E>1 kcV) X-rays can penetrate a large gas column 
(Nh > 10 22 cm" 2 ) forming an X-ray Dissociation Region 
(XDR: iMalonev et "all ([1996)). The key quantity for un- 
derstanding the XDR heating and cooling is the ratio 
of the X-ray energy deposition rate per particle (set by 
the X-ray luminosity and attenuating column) to the gas 
density which sets the recombination and cooling rates. 
We do not have a direct view of the AGN in the Clover- 
leaf. The observed X -ray emission is believed to be scat- 
tered ( Chartas 2000) , and even this scattered emission is 
subject to an uncertain absorption in the gas along the 
line of sight. Thus the true X-ray flux must be inferred 
or modeled based on the observed spectrum. At the 
bolometric luminos ity estimated for the Cloverleaf AGN 
(Lboi ~ 7 x 10 13 L^(lLutz et al.l l200lK the 1- 20 keV X- 



rays typically carry ~ 5% of £h n i ()Mushotzkv et al.lll993f) 
(though with a scatter of a factor o f ~ 3 in this relation) . 
This implies that L x « 3 x 1O 12 L . iChartas et all (|2007h 
present a model (their Figure 9) explaining the Clover- 
leaf X-ray spectra from Chandra and XMM-Newton in 
which the total 2-10 keV luminosity of the central region 
is 2 x 10 46 ergs _1 or 5x 10 12 L Q (instrinsic, with m ~ 10), 
based on an observed flux which is a factor of ^200 lower. 
While these inferences clearly have large uncertain- 
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Fig. 6. — CO abundance (relative to H nuclei) in the distance- 
gas density space of the XDR model. At small distances and low 
densities, gas is largely atomic. The blue contours show gas tem- 
perature. 

ties, it is very likely that the hard X-ray luminosity ap- 
proaches the 40-120/im far-infrared luminosity, and 
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The large column densities derived in the CO-line anal- 
ysis mean that much of this hard X-ray emission will be 
absorbed. Unlike in a PDR, the X-rays input a much 
greater fraction their energy (~0.1-1) into the gas than 
the dust, and since the CO lines will be an important 
coolant, this ratio immediately indicates that the AGN 
could readily power the CO emission we observe. 

To explore this possibility, we have generated XDR 
models for parameters appropriate to the Cloverleaf sys- 
tem. The input parameters are the total gas density and 
the incident X-ray flux, which is set (for fixed Lx) by the 
distance to the AGN and the column density of absorb- 
ing gas between the hard X-ray source and the modeled 
region. The results also depend on the optical depth in 
the coolings lines; we have used a total hydrogen column 
per unit linewidth of 10 21 cm -2 /kms -1 , as suggested by 
the CO line analysis. 

Figure [5] shows the predicted CO abundance as a func- 
tion of gas density and cloud distance (denoted R) from 
the nucleus. The assumed X-ray luminosity is conserva- 
tively adopted to be 10 46 ergs -1 and the total attenuat- 
ing column is Nn, a tt = 3 x 10 23 cm -2 , of order but less 
than the estimated total column from the CO analysis. 
This results in a factor of 22 attenuation and the result- 
ing X-ray flux at 600 pc distance is thus 10 ergs -1 cm -2 . 
Larger values of Lx or smaller values of N at t shift the con- 
tours up and to the right. We have adopted solar abun- 
dances for elements in the gas-phase. At small R and 
relatively low density (lower left corner of the plot), the 
CO abundance is small, as the gas is warm and atomic. 
However, over most of the plotted range, the CO abun- 
dance is large (more than ~a few x 10 -5 ), reaching nearly 
3 x 10 -4 at the largest R and densities, for which all the 
gas-phase carbon is in CO. 

As noted above, the surface brightness of the CO emis- 
sion in the Cloverleaf is very large. The estimated total 
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Fig. 7. — CO rotational cooling for X-ray heated clouds in the 
XDR model space. The contours show surface brightness of all the 
CO lines in units of erg / s / cm 2 . Values of 0.2 and 0.6 correspond 
to the observed surface brightness assuming filling factors between 
1/3 and 1 in the VS03 disk. Dashed blue contours show thermal 
pressure, values denoted with lognT 

CO rotational line luminosity (equation [2]), distributed 
over the VS03 disk with unit area filling factor implies a 
surface flux of approximately 0.2 ergs cm -2 s -1 . If the 
area filing factor is less than unity, then even larger local 
surface brightness is required. We plot this local CO sur- 
face brightness predicted by our XDR models in Figure[71 
and conclude that for densities above 3 x 10 4 cm -3 , and 
R from 500-1500 pc, our model is consistent with the 
observed CO cooling, and suggests area filling factors in 
the range of 0.3-1. Of course, even if the area filling fac- 
tor is unity, we note that the volume filling factors can 
still be quite s mall. Our result s are c onsistent with the 
XDR model of iMeiierink et alj (|2007l ). who present the 
CO J = 7 ^ 6 to J = 3 ^ 2 intensity ratio as a function 
on incident X-ray flux and density (their Figure 6). The 
observed ratio is 8.9-12.3, depending whichJ = 3^2 
measurement is adopted, suggesting an X-ray flux of 5- 
20 ergs -1 cm -2 , which compares well with our working 
value of ~10 at d = 600 pc. 

More generally, the XDR models indicate that CO ro- 
tational line emission can be ~ 0.1 of the absorbed X-ray 
energy, so only ~ 1% of the hard X-ray luminosity of the 
AGN needs to be absorbed in a CO-luminous XDR to 
produce the CO emission. However, we emphasize that 
while the XDR can easily heat the gas, it is not likely 
to appreciably heat the co-extant dust the way that the 
PDR front would. In an XDR, the ratio of total CO 
rotational line flux to the locally- generated far-infrared 
flux from grains (due to absorption of emitted line pho- 
tons and to UV photons produced by excitation and re- 
combination) can easily exceed 10 -3 , and can be as large 
as ~ 0.1. (This ratio depends on the cloud column den- 
sity as well as the column density per unit linewidth, 
since the far-infrared dust emission depends on the for- 
mer.) Since the observed ratio is ~ 5 x 10 -4 , the observed 
far-IR / submillimeter dust emission can not be pow- 
ered by the XDR processes. There must be additional 
source(s) of local grain heating, e.g., star formation, or 
re-radiated continuum from the AGN. 
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Fig. 8.— Ratio of the total CO cooling to the [CI] J = 1 —* 
cooling in the XDR model. The heavy contour corresponds to the 
observed value (130). 

As with the PDR, the XDR framework does not per- 
fectly account for the observed [CI] intensities, though 
it may provide a better match in terms of total cooling. 
The XDR model generally predicts lower [CI] relative to 
CO than the PDR, and the [CI] J = 1 -> line is well- 
matched with the XDR model. The ratio of the total 
CO emergent intensity including all transitions to the 
[CI] J = 1 — > intensity is plotted in Figure [51 The ob- 
served value of 130 is highlighted, as are ratios a factor of 
two above and below this value to allow for systematic 
uncertainties in both. As with the analysis described 
above, the allowed region corresponds to sub-kpc dis- 
tances and densities greater than 10 5 cm~ 3 . However, it 
remains true that the [CI] temperature is far too low; the 
[CI] J = 2— >1 to J=l— >0 ratio (R = 2.3 in power units 
per W03, W05) suggests a physical temperature of only 
~30 K, a value generally inconsistent with both the the 
XDR and the PDR models. W05 suggest that the excess 
[CI] J=l — >0 emission may arise in an (additional) cool 
gas component exterior to the material traced with the 
CO J = 7 — > 6. Such a scenario may be consistent with 
our findings, as we allow for a second component of low- 
excitation material outside the disk. However, this would 
mean that the [CI] J = 1 — > which is associated with 
the mid-J CO would be less, corresponding to higher ra- 
tios in Figure [SJ requiring higher densities and lower dis- 
tances from the central source. A large allocation of the 
[CI] J = 1 — > to this external component is difficult to 
reconcile with the overall XDR picture. Another possi- 
bility is that the measurements of the relatively weak [CI] 
intensities are subject to systematic errors. We do note 
that our measurement favors a slightly higher flux for the 
[CI] J = 2^ 1 than W05 (giving R ~ 3.8, for T-40 K), 
though with low significance. We do note that the W05 
measurement of [CI] jtwo puts the [CI] ratio well below 
that of M82 (R = 4.3, per fStutzki et alj|1997[ )). the only 
external galaxy for which the J = 2 — * 1 measurement is 
published. 

In the context of the XDR model, the above analyses 
suggest that the gas density is high: nn 2 > 5 x 10 4 cm~ 3 , 
for the bulk of the material. This density is higher 



than required by the CO intensity ratios alone (Fig- 
ure [3]); it derives in essence from requiring the total 
observed CO luminosity to be generated in a volume 
constrained by the VS03 disk model. High densities 
have been suggested given the measurement of high- 
dipole-moment s pecie s, beginning with HCN J = 1 — * 
Solomon etal] (|2003[) '). and now with HCO+ J = 1 -> 
Riechers et al.ll2006a[) . The two transitions are compa- 
rable: L HCO + = 0.8±0.2Lhcn, arid the fractional inten- 
sity relati ve to Lpm nicely match es the linear correlation 
found by iGao fc Solomon! (2004) for the local-universe 
LIRGS and ULIRGS. While mid-IR pu mping schemes 
could p otentially account for these lines, iRiechers et al.l 
(2006a) conclude that the match between the HCN and 
HCO + lines is best explained by both species produc- 
ing optically-thick thermalized emission (at least for the 
J=1^0 transition), requiring njj 2 > 10 5 . This is fully 
consistent with the both the XDR and the high-density 
PDR scenarios. 

4.4.3. Turbulence, Cosmic Rays 

Finally, we note that there are other potential energy 
sources which will preferentially heat the gas relative to 
the dust. Cosmic-ray-ionization heating has been pro- 
posed to explain the warm m olecular material in the 
nucle i of M 82 and NGC 2 53 (iHailev-Dunsheath et all 
[20081: Bradford et all l200l ISuchkov et al l 11993) . Cos- 
mic Rays are generated in supernovae associated with 
massive star formation, and like X-rays, they penetrate 
throughout the bulk of the molecular material. Sim- 
ple dissipation of me chanical energy is als o a poten- 
tial heating source. Bradford et all (120051) apply the 
shock models oflDraine et al.l (I1983D : iDraine fc Roberge! 
(|1984D : iRoberge fc Draind (| 19901) to the molecular gas 
around SGR A* and find that low- velocity (v ~10- 
20 kms -1 ) magnetohydrodynamic (MHD) shocks pro- 
duced by cloud-cloud collisions can produce warm ma- 
terial which cools primarily via the mid-J CO rotational 
transition, again without appreciable heating of the dust. 
The arguments presented in these studies would apply to 
the Cloverleaf host, and could also boost the mid-J CO 
emission relative to that which arises in the PDR alone. 

4.4.4. Prospects for Atomic Gas Tracers 

Further distinction between the PDR scenario and the 
X-ray and other various bulk heating scenarios could be 
provided by the atomic fine-structure transitions. In par- 
ticular, the [OI] 63/im and [CII] 158/im will be impor- 
tant. The ratio of these fine-structure line intensities to 
the total CO cooling is a powerful discriminant: in the 
PDR the fine structure lines carry an order of magnitude 
more power than the CO since most of the UV heating 
is deposited in the atomic gas, while in the XDR, the 
energy is more uniformly distributed between the atomic 
and molecular components. Indeed, our models suggest 
that the [OI] 63/im intensity and the total CO rotational 
cooling are comparable for the XDR. With the Clover- 
leaf's unusual luminosity and lensing factor, these fine- 
structure lines may be accessible with the spectrometers 
on Herschel. 

4.4.5. Conclusion - A Composite Solution 

We conclude that there are several potential gas heat- 
ing sources, each of which could contribute substantially 
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to powering the CO emission in the Cloverleaf host. The 
far-IR / submillimeter continuum component appears 
distinct from the mid-IR component, and the PAH frac- 
tional luminosities relative to this far-IR dust component 
are similar to those observed in starburst systems. These 
facts are difficult to explain without a massive starburst. 
On the other hand, since the fractional energy carried 
in the molecular gas exceeds that of the local starburst 
analogs by factors of 2-5, and there is a powerful X-ray 
source which is demonstrably capable of contributing 10- 
100% of the required heating, it is difficult to exclude the 
X-rays as a heating source. We conclude that the molec- 
ular material around the Cloverleaf QSO is most likely 
heated by both UV photons from young stars and and 
X-rays from the accretion zone, with neither dominant 
by more than ~1 order of magnitude. 

X-rays which penetrate and heat the bulk of the molec- 
ular gas will impact the properties of the star formation, 
particularly the stellar mass function. Simple thermo- 
dynamic arguments along the lines of the Jeans mass 
or Bonnor-Ebert mass show that the characteristic mass 
required for gravitational collapse increases as the 'mini- 
mum temperature' of the gas increases. This is precisely 
the effect that a bulk heating mechanism would produce. 
Neglecting the effects of magnetic fields and rotation, this 
mass is given by 

M be = pl/2 G3/2 oc CX pj^, (5) 

where T, n, and P are temperature, number density, 
and thermal pressure, respectively of the environment 
from which the stars must form, and Cbe is a numer- 
ical constant. Detailed theoretical approaches confirm 
the sense of this relation, finding a scaling between the 
characteristic mass scale M* in a stellar IMF, and the 
minimum temperature to which molecular gas can cool — 
M* cx T^ in — with exponen t 7 ranging from 1.7 (obtained 
in numerical experiments (jJappsen et al J 120051)1 t o 3.35 
(obtained in an analytic treatment (|Larsonl 1198511. In 
the Galaxy, M* ~ 0.5 M and T min ~ 8 K, values nu- 
merically consistent with Equation O While our anal- 
ysis does not directly measure T m j n , it must be true 
that T m i n is increased relative to the Galaxy. A plau- 
sible assertion is that T min scales as the typical tem- 
perature derived from fits to CO line ratios: at least 
~50K per our analysis of the Cloverleaf fluxes, compared 
with -22 K in the inner Galaxy p er the COBE FIRAS 
measurements (|Fixsen et al.l Il999f) . This simple scaling 
would suggest that M* in the Cloverleaf starbur st is 4- 
15 tim es the Galactic value: some 2-5 M Q . As iLarsor] 
(1998) and others have pointed out, such a top-heavy 
IMF converts a given mass of gas into a greater total stel- 
lar (and thus far-IR) luminosity than with the Salpeter 
IMF. This is the leading scenario proposed to explain the 
factor of ~3-5 discrepancy between the observed stel- 
lar mass buildu p and the star-formation history in the 
4 < z < 1 era (iPerez-Gonzalez et aLll2008t iDavel 120081 : 
iHopkins fe; Beacomll2006l ). While the Cloverleaf may be 
somewhat more extreme than the typical star-forming 
system, the increased prevalence of AGN in this early 
epoch may generally result in more massive stars than 
those formed in the Galaxy today. 

4.5. Additional Features in the Spectrum 



While it has yet to be explored completely in a extra- 
galactic source, the short submillimeter band is expected 
to host low-lying transitions of light molecules and ions 
other than CO, some of which may be bright. We now 
discuss some tentative identifications which could be fol- 
lowed up with the large collecting area and spectral res- 
olution of the interferometers. 

4.5.1. Water in the Cloverleaf 

The feature at 277.6 GHz is fit with a Gaussian line of 
^rcst = 987.4±0.55GHz, adopting the systemic redshift of 
the Cloverleaf. This might be identified with the lowest- 
lying transition of ortho-hydronium o — HsO + (Op — > 1q) 
t'rcst = 984 GHz), but unlike the other low-lying HsO + 
transitions (v les t — 396, 388, 364 GHz), it is not expected 
to be bright given the lack o f radiative pumping pathways 
to excite the upper level ([Phillips et al.l fl992f ). More- 
over, the abundance of H3O 4 " even in extreme Arp 220- 
like nuclei is believed to be less than 1/10 that of wa- 
ter ([van der Tak et al.l l20081. A better match to the fit, 
and a more convincing scenario is that this feature is the 
987.9 GHz (2q,2 — ► li,i) transition connecting the first 
excited level of p-H 2 at 54 K to the second excited 
level at 101 K. Attempts to detect hig h-redshift rota- 
tional water transitions ha ve been made. lEncrenaz et al.l 
(19931: [Casoli et ail (|1994D report tentative detections of 
the 2i,i — > 2 0i2 transition (u Test = 752.0 GHz) in IRAS 
F10214+4724 (z = 2.3). Remarkably, the reported in- 
tensity (0.65 Kkms -1 at 30-m) is only a factor of 0.45 
times that of the nearb y CO J = 6 — > 5 t ransition 
(1.4 Kk ms" 1 30-m. per ISolomon et all (|1992f) 1. More 
recently, (|Riechers et al.ll2006bt ) report an upper limit on 
the 3i, 3 -> 2 2 , (183 GHz) transition in MG 0751+2716 
at z = 3.2. 

The a priori interpretation of water spectra is diffi- 
cult even with multiple line detections, as evidenced by 
the work analyzing the ISO spectra toward SGR B2 
(jNeufeld et a ll 119951: ICernicharo eta l. 2006 b|) and the 
Orion outflow (|Harwit et alj Il998t ICernicharo et al.l 
2006a). Morever, since most of the low-lying water lines 
are not accessible at zero redshift from the ground, there 
are no good Galactic templates with which to compare 
the 20,2 - * I1.1 measurement. Surveys have been con- 
ducted only in the ground-state li — > In 1 (538 p,m) 
line w it h the KAO, S WAS, a nd ODIN (e . g. lAlihbv et all 
l|200lft : lNeufeld et all (l2003h : ISnell etal] (boOOD I. finding 
low abundances (< 10~ 8 1 or upper limits in GMCs and 
cloud cores. Similarly, IWilson et al.l (|2007f ) derive up- 
per limits to the water abundance of < 10~ 8 from non- 
detections in nearby starburst galaxies, albeit with the 
large (2.1') ODIN beam. However, much larger water 
abundances (more than 10 -5 ) are inferred in warm cores 
and outflows, notably the Orion outflow. 

In spite of the ODIN ground-state non-detections 
in local starburs ts, the ISO far-IR spectrum of the 
Arp 220 nucleus ([Gonzalez- Alfonso et al] [20041 ) demon- 
strates that water can be abundant and produce pow- 
erful features even in a large-beam (i.e. average) ex- 
tragalactic spectrum, at least in extreme sources. The 
179 /im absorption from the nucleus shows an equivalent 
width of ~400kms _1 (in absorption). These investiga- 
tors quote H 2 column densities of 2-10x 10 17 cm~ 2 to- 
ward the Arp 220 nucleus, corresponding to abundances 
of l-5xl0 -8 for Nh2 ~ 10 25 cm~ 2 . We measure an 
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equivalent width of ~580kms , which is of the same 
order, albeit in emission. We do note that the models of 
ICernicharo et alj (|2006aj ) suggest that even in the SGR 
B2 geometry, the 2o,2 - * li.i transition in particular could 
be brought into emission (in contrast with most of the 
other transitions) due to a favorable pumping / cascade 
network in the para levels. While there is no straight- 
forward detailed reconciliation of the tentative detection 
with published models, we note that water emission is 
more easy to understand in the context of an XDR than 
in a traditional PDR. Our XDR model suggests a water 
abundance of order 10~ 7 , greater than that derived in 
Arp 220, and more than would be expected in the bulk 
of molecular material of a PDR. 

4.5.2. Upper limits, absorption identification 

Before concluding, we briefly note for completeness 
upper limits for two other transitions, marked in the 
spectrum in Figure [T] The radical CH + has its ground 
state transition ( J = 1 — * 0) at z/ rcst = 806GHz (z^ bs = 
234.7GHz). The spectrum shows this channel is above 
the continuum at the 1.5a level; a formal 2.5a upper 
limit to the flux is 14 Jy km/s. The LiH J = 2^1 transi- 
tion at fc'rest = 887GHz (fobs = 249.4GHz) is also above 
the continuum at low significance, with a formal upper 
limit of 12.3 Jy kms -1 . 

The absorption feature at 287 GHz is significant, 
at 3cr, Interestingly, this frequency corresponds to the 
557 GHz ground state of 0-H2O at z=0.93, perfectly 
within the redshift r ange expected for the Cloverleaf lens 
(jKneib et alj Il998bl ) (z = 0.9 ± 0.1). The equivalent 
width of the absorption is ^600 km/s, which would im- 
ply that there must be nearly complete absorption over a 
several hundred km/s velocity range. A galaxy perfectly 
positioned along the line of sight could potentially do 
this. For r to be unity over 600 km/s requires a molecular 
hydrogen column of 6 x 10 23 cm~ 2 , assuming an ortho-to- 
para ratio of unity, a water abundance of 10 -8 , and that 
there is negligible population in the upper levels. The 
lens is located 0.6 arcsec from the projected Cloverleaf 
images, thus the light passes within <5 kpc of the lens, 
so the geometry is not implausible. Further investigation 
of this intriguing possibility will require higher spectral 
and spatial resolution measurements. 

5. CONCLUSIONS 

We have built and are now using the first broadband 
spectrograph for the millimeter / submillimcter, provid- 
ing a new observational approach for studying galaxies 
at all redshifts. Among the first experiments is the mea- 
surements of the rest-frame 272-444 fim spectrum of the 
Cloverleaf system which contains the bolometric peak in 
the CO rotational spectrum. In our analysis of the CO 
intensities and other published data we find the follow- 
ing: 

1. The intensity ratios of the J = 3 — > 2 through 
J = 9 — > 8 CO transitions can be produced with a 
single gas component. The excitation is high, with 
thermal pressure > 10 6 KcitT 3 . The 0.2-5xl0 10 
M© of excited molecular gas in the Cloverleaf host 
must fi ll a projected size of at least half the 650-pc- 
radius iVenturini fc Solomon! (|2003l ) disk. A more 



compact distribution is not physical based on the 
ab solute line intensit ies . 
2. As lLutz et al.1 (|2007T ) conclude, the match between 
the far-IR dust component and the PAH emission 
suggests that far-IR-emitting dust component rep- 
resents a massive starburst. However, the molecu- 
lar gas cooling in the CO transitions alone is a large 
fraction (6 x 10 -4 ) of the bolometric luminosity of 
this starburst. This is a factor of a few larger than 
in the local-universe starburst galaxies for which 
mid-J CO intensities are available. 



3. Given the powerful AGN of the Cloverleaf system, 
X-rays are likely to be an important energy source 
for molecular gas even several hundred parsecs from 
the nucleus. Our XDR model show that ~5% of 
the bolometric luminosity of the AGN in hard X- 
rays can easily provide the energy input required 
to match the observed CO cooling. The X-rays 
do not appreciably contribute to the far-IR dust 
continuum. 

Our interpretation is that the Cloverleaf host is indeed 
a undergoing a massive starburst, but that it has addi- 
tional energy input into the ISM via hard X-rays orig- 
inating in the accretion zone. As a result, the stellar 
IMF is likely biased toward much higher masses than 
in the Galaxy, with a typical stellar mass M* of sev- 
eral M©. With this characteristic mass replacing the 
Salpeter value of ~0.5 M©, the gas consumption rate 
at the observed luminosity is several t i mes l ower than 
that given by applying the iKennicuttl (|1998| ) prescrip- 
tion. The starburst we are witnessing in the Cloverleaf 
may ther efore extend much longer than the 30 Myr esti- 
mated bv lLutz et a l. (2007); it could last for a few hun- 
dred Myr. A similar scenario, though perhaps less ex- 
treme, may be typical of the star-forming galaxies in the 
first half of the Universe's history when energy release 
peaked. 



We are indebted to the staff of the Caltech Submil- 
limeter Observatory for their help in Z-Spec's commis- 
sioning and observing. We acknowledge Peter Ade and 
his group for some of the filters and Lionel Duband for 
the 3 He / 4 He refrigerator in Z-Spec, and are grateful for 
their help in the early integration of the instrument. We 
benefitted from conversations with Tom Phillips, Paul 
Goldsmith, Simon Radford, and Andy Harris, as well as 
helpful comments from Xinyu Dai and an anonymous 
referee. Finally, we acknowledge the following grants 
and fellowships: NASA SARA grants NAGS-11911 and 
NAGS-12788, an NSF Career grant (AST-0239270) and 
a Research Corporation Award (RI0928) to J. Glenn, 
a Caltech Millikan and JPL Director's fellowships to 
C.M.B., a NRAO Jansky fellowship to J. Aguirre, and 
NASA GSRP fellowship to L. Earle. The research de- 
scribed in this paper carried out at the Jet Propulsion 
Laboratory, California Institute of Technology, was done 
so under a contract with the National Aeronautics and 
Space Administration. 



12 



Bradford et al. 



REFERENCES 



Alloin, D., Guilloteau, S., Barvainis, R., Antonucci, R., & Tacconi, 
L. 1997, A&A, 321, 24 

Ashby, M. L. N., Bcrgin, E. A., Plume, R., Carpenter, J. M., 
Neufeld, D. A., Chin, G., Erickson, N. R., Goldsmith, P. F., 
Harwit, M., Howe, J. E., Kleiner, S. C, Koch, D. G., Patten, 
B. M., Schieder, R., Snell, R. L., Stauffer, J. R., Tolls, V., Wang, 
Z., Winnewisser, G., Zhang, Y. F., & Melnick, G. J. 2000, ApJ, 
539, L115 

Barvainis, R., Antonucci, R., & Coleman, P. 1992, ApJ, 399, L19 
Barvainis, R., Antonucci, R., Hurt, T., Coleman, P., & Reuter, 

H.-P. 1995, ApJ, 451, L9+ 
Barvainis, R., Maloney, P., Antonucci, R., & Alloin, D. 1997, ApJ, 

484, 695 

Barvainis, R., Tacconi, L., Antonucci, R., Alloin, D., & Coleman, 
P. 1994, Nature, 371, 586 

Bradford, C. M., Ade, P. A. R., Aguirre, J. E., Bock, J. J., 
Dragovan, M., Duband, L., Earle, L., Glenn, J., Matsuhara, 
H., Naylor, B. J., Nguyen, H. T., Yun, M., & Zmuidzinas, J. 
2004, in Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, Vol. 5498, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, ed. C. M. 
Bradford, P. A. R. Ade, J. E. Aguirre, J. J. Bock, M. Dragovan, 
L. Duband, L. Earle, J. Glenn, H. Matsuhara, B. J. Naylor, H. T. 
Nguyen, M. Yun, & J. Zmuidzinas, 257-+ 

Bradford, C. M., Nikola, T., Stacey, G. J., Bolatto, A. D., Jackson, 
J. M., Savage, M. L., Davidson, J. A., & Higdon, S. J. 2003, ApJ, 
586, 891 

Bradford, C. M., Stacey, G. J., Nikola, T., Bolatto, A. D., Jackson, 

J. M., Savage, M. L., & Davidson, J. A. 2005, ApJ, 623, 866 
Brown, R. L., & Vanden Bout, P. A. 1992, ApJ, 397, L19 
Casoli, F., Gerin, M., Encrenaz, P. J., & Combes, F. 1994, A&A, 
287, 716 

Cernicharo, J., Goicoechea, J. R., Daniel, F., Lerate, M. R., Barlow, 
M. J., Swinyard, B. M., van Dishoeck, E. F., Lim, T. L., Viti, 
S., & Yates, J. 2006a, ApJ, 649, L33 

Cernicharo, J., Goicoechea, J. R., Pardo, J. R., & Asensio-Ramos, 
A. 2006b, ApJ, 642, 940 

Chartas, G. 2000, ApJ, 531, 81 

Chartas, G., Eracleous, M., Agol, E., & Gallagher, S. C. 2004, ApJ, 
606, 78 

Chartas, G., Eracleous, M., Dai, X., Agol, E., & Gallagher, S. 2007, 

ApJ, 661, 678 
Dave, R. 2008, MNRAS, 385, 147 

Draine, B. T., & Roberge, W. G. 1984, ApJ, 282, 491 
Draine, B. T., Roberge, W. G., & Dalgarno, A. 1983, ApJ, 264, 
485 

Earle, L., Ade, P., Aguirre, J., Aikin, R., Battle, J., Bock, J., 
Bradford, C. M., Dragovan, M., Duband, L., Glenn, J., Griffin, 

G. , Hristov, V., Maloney, P., Matsuhara, H., Naylor, B., Nguyen, 

H. , Yun, M., & Zmuidzinas, J. 2006, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, 
Vol. 6275, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series 

Encrenaz, P. J., Combes, F., Casoli, F., Gerin, M., Pagani, L., 

Horellou, C, & Gac, C. 1993, A&A, 273, L19+ 
Fixsen, D. J., Bennett, C. L., & Mather, J. C. 1999, ApJ, 526, 207 
Frayer, D. T., Koda, J., Pope, A., Huynh, M. T., Chary, R.- 
R., Scott, D., Dickinson, M., Bock, D. C.-J., Carpenter, J. M., 
Hawkins, D., Hodges, M., Lamb, J. W., Plambeck, R. L., Pound, 
M. W., Scott, S. L., Scoville, N. Z., & Woody, D. P. 2008, ApJ, 
680, L21 

Gao, Y., & Solomon, P. M. 2004, ApJ, 606, 271 
Gonzalez-Alfonso, E., Smith, H. A., Fischer, J., & Cernicharo, J. 

2004, ApJ, 613, 247 
Giisten, R., Philipp, S. D., Weifi, A., & Klein, B. 2006, A&A, 454, 

L115 

Hailey-Dunsheath, S., Nikola, T., Stacey, G. J., Oberst, T. E., 

Parshley, S. C, Bradford, C. M., Ade, P. A. R., & Tucker, C. E. 

2008, ArXiv e-prints 
Harwit, M., Neufeld, D. A., Melnick, G. J., & Kaufman, M. J. 1998, 

ApJ, 497, L105+ 
Hazard, C, Morton, D. C, Terlevich, R., & McMahon, R. 1984, 

ApJ, 282, 33 

Henkel, C, Mauersberger, R., Wiklind, T., Huettemcistcr, S., 

Lcmmc, C, & Millar, T. J. 1993, A&A, 268, L17 
Hopkins, A. M., & Beacom, J. F. 2006, ApJ, 651, 142 



Inami, H., Bradford, M., Aguirre, J., Earle, L., Naylor, 
B., Matsuhara, H., Glenn, J., Nguyen, H., Bock, J. J., 
Zmuidzinas, J., & Ohyama, Y. 2008, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, 
Vol. 7020, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series 

Jappscn, A.-K., Klessen, R. S., Larson, R. B., Li, Y., & Mac Low, 
M.-M. 2005, A&A, 435, 611 

Kaufman, M. J., Wolfire, M. G., & Hollenbach, D. J. 2006, ApJ, 
644, 283 

Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 

Kneib, J.-P., Alloin, D., Mellier, Y., Guilloteau, S., Barvainis, R., 

& Antonucci, R. 1998a, A&A, 329, 827 
Kneib, J.-P., Alloin, D., & Pello, R. 1998b, A&A, 339, L65 
Larson, R. B. 1985, MNRAS, 214, 379 
— . 1998, MNRAS, 301, 569 

Lutz, D., Sturm, E., Tacconi, L. J., Valiante, E., Schweitzer, M., 
Nctzcr, H., Maiolino, R., Andrcani, P., Shcmmer, O., & Vcillcux, 
S. 2007, ApJ, 661, L25 

Maloney, P. R., Hollenbach, D. J., & Tielens, A. G. G. M. 1996, 
ApJ, 466, 561 

Meijerink, R., Spaans, M., & Israel, F. P. 2007, A&A, 461, 793 
Mushotzky, R. F., Done, C, & Pounds, K. A. 1993, ARA&A, 31, 
717 

Naylor, B. J., Ade, P. A. R., Bock, J. J., Bradford, C. M., Dragovan, 
M., Duband, L., Earle, L., Glenn, J., Matsuhara, H., Nguyen, 
H., Yun, M., & Zmuidzinas, J. 2003, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Series, 
Vol. 4855, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. T. G. Phillips & J. Zmuidzinas, 
239-248 

Neufeld, D. A., Bcrgin, E. A., Melnick, G. J., & Goldsmith, P. F. 

2003, ApJ, 590, 882 
Neufeld, D. A., Lcpp, S., & Melnick, G. J. 1995, ApJS, 100, 132 
Papadopoulos, P. P., Isaak, K. G., & van der Werf, P. P. 2007, ApJ, 

668, 815 

Perez-Gonzalez, P. G., Rickc, G. H., Villar, V., Barro, G., Blaylock, 
M., Egami, E., Gallego, J., Gil de Paz, A., Pascual, S., Zamorano, 
J., & Donley, J. L. 2008, ApJ, 675, 234 

Phillips, T. G., van Dishoeck, E. F., & Keene, J. 1992, ApJ, 399, 
533 

Riechers, D. A., Walter, F., Carilli, C. L., Weiss, A., Bertoldi, F., 
Menten, K. M., Knudsen, K. K., & Cox, P. 2006a, ApJ, 645, L13 

Riechers, D. A., Weiss, A., Walter, F., Carilli, C. L., & Knudsen, 
K. K. 2006b, ApJ, 649, 635 

Roberge, W. G., & Draine, B. T. 1990, ApJ, 350, 700 

Snell, R. L., Howe, J. E., Ashby, M. L. N., Bergin, E. A., Chin, 
G., Erickson, N. R., Goldsmith, P. F., Harwit, M., Kleiner, S. C, 
Koch, D. G., Neufeld, D. A., Patten, B. M., Plume, R., Schieder, 
R., Stauffer, J. R., Tolls, V., Wang, Z., Winnewisser, G., Zhang, 
Y. F., & Melnick, G. J. 2000, ApJ, 539, L101 

Solomon, P., Vanden Bout, P., Carilli, C, & Guclin, M. 2003, 
Nature, 426, 636 

Solomon, P. M., Downes, D., & Radford, S. J. E. 1992, ApJ, 398, 
L29 

Solomon, P. M., & Vanden Bout, P. A. 2005, ARA&A, 43, 677 
Stacey, G. J., Hailey-Dunsheath, S., Nikola, T., Oberst, T. E., 
Parshley, S. C, Benford, D. J., Staguhn, J. G., Moseley, 
S. H., & Tucker, C. 2007, in Astronomical Society of the 
Pacific Conference Series, Vol. 375, From Z-Machincs to ALMA: 
(Sub)Millimeter Spectroscopy of Galaxies, ed. A. J. Baker, 
J. Glenn, A. I. Harris, J. G. Mangum, & M. S. Yun, 52-+ 
Stutzki, J., Graf, U. U., Haas, S., Honingh, C. E., Hottgcnroth, D., 
Jacobs, K., Schieder, R., Simon, R., Staguhn, J., Winnewisser, 
G., Martin, R. N., Peters, W. L., & McMullin, J. P. 1997, ApJ, 
477, L33+ 

Suchkov, A., Allen, R. J., & Heckman, T. M. 1993, ApJ, 413, 542 
van der Tak, F. F. S., Aalto, S., & Meijerink, R. 2008, A&A, 477, 
L5 

van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & 

van Dishoeck, E. F. 2007, A&A, 468, 627 
Venturini, S., & Solomon, P. M. 2003, ApJ, 590, 740 
Ward, J. S., Zmuidzinas, J., Harris, A. I., & Isaak, K. G. 2003, 

ApJ, 587, 171 

Weifi, A., Downes, D., Henkel, C, & Walter, F. 2005, A&A, 429, 
L25 



Warm Molecular Gas in the Cloverleaf 



13 



WeiB, A., Hcnkcl, C, Downes, D., & Walter, F. 2003, A&A, 409, Yun, M. S., Scoville, N. Z., Carrasco, J. J., & Blandford, R. D. 

L41 1997, ApJ, 479, L9+ 

Wilner, D. J., Zhao, J.-H., & Ho, P. T. P. 1995, ApJ, 453, L91+ 
Wilson, C. D., Booth, R. S., Olofsson, A. O. H., Olberg, M., 

Pcrsson, C. M., Sandqvist, A., Hjalmarson, A., Buat, V., 

Encrenaz, P. J., Fich, M., Frisk, U., Gerin, M., Rydback, G., 

& Wiklind, T. 2007, A&A, 469, 121 



